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Abstract 

We present a simple and accessible method which uses contour integration 
methods to derive formulae for functional determinants. To make the presen- 
tation as clear as possible, the general idea is first illustrated on the simplest 
case: a second order differential operator with Dirichlet boundary conditions. 
The method is applicable to more general situations, and we discuss the way in 
which the formalism has to be developed to cover these cases. In particular, we 
also show that simple and elegant formulae exist for the physically important 
case of determinants where zero modes exist, but have been excluded. 



1 Introduction 



The need to calculate functional determinants, usually of second order differential op- 
erators, arises in many different fields. For example, the leading order contribution 
to path-integrals often involves the integration of the exponential of a quadratic form, 
which can be carried out exactly, and yields a functional determinant |2] . The ex- 
plicit formulae for one-dimensional functional determinants of this kind are not simply 
of interest because of their utility, but also because rather general results may be ob- 
tained. For instance, the functional determinant of an operator L can often be given 
explicitly in terms of the solutions of the homogeneous equation Ly = 0. The elegance 
and generality of the results that can be obtained has attracted mathematicians, es- 
pecially in the last decade or so. However, many researchers who are simply searching 
for an explicit form for a particular differential operator may find these treatments 
rather abstract, and, in addition, some will be interested in the value of these deter- 
minants with the zero modes extracted. These zero modes will usually exist whenever 
a quadratic form in a path-integral is the result of performing a linearisation about 
a solution with non-trivial spatial or temporal dependence. By Goldstone's theorem, 
this symmetry breaking will give rise to zero modes, which must be excluded from the 
quadratic form if this is to be a consistent leading order approximation [3 j. Thus the 
functional determinant obtained by evaluating the path-integral is as before, but with 
these modes absent. 

This paper is designed to give a simple, and easily accessible, derivation of the 
main formulae for functional determinants, and also to show how these results may be 
modified to give analogous expressions for determinants with the zero mode (or modes) 
extracted. We will avoid unnecessary mathematical abstraction, concentrate on the 
most straightforward results and give examples to illustrate the methods; a forthcoming 
paper will discuss some of the technical points in more detail jlj. Specifically, we will 
use only contour integration methods and elementary facts concerning the solutions of 
differential equations and concentrate on second order differential operators. An earlier 
study of determinants with the zero modes extracted |Hj obtained results by invoking 
a regular isat ion procedure in which the zero modes became non-zero, extracting them, 
and then removing the regularisation. The problem with this approach is that it is 
not clear that the result is independent of the regularisation. Within the contour 
integration adopted here, the zero modes are not included from the beginning, and 
there is no need to introduce any regularisation procedure for them. 

An early review of functional integrals and their evaluation [6 3 gives a prescription 
for finding the functional determinant of a second order differential operator with no 
first order derivative and with Dirichlet boundary conditions. The method consists of 
discretising the functional integral, explicitly carrying out the resulting multidimen- 
sional Gaussian integral, and showing that, on return to the continuum, the determi- 
nant satisfies the homogeneous equation with given boundary conditions. This method 
is also clearly presented in Schulman's book on path-integrals |2J. During the 1970's 
quantum field theory in general, and functional integral techniques in particular, be- 
came very popular and many studies of non-trivial field configurations such as solitons 
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and instantons appeared. In some cases [7] the calculation of the determinants in- 
volved in these studies were carried out by explicitly determining the eigenvalues of 
the Sturm-Liouville problem and calculating their product. This had the advantage 
that the omission of the zero mode required in these cases could be carried out explic- 
itly. On the other hand, Coleman, in an appendix to his 1977 lectures on instantons at 
Erice jS], gives the outline of a rather straightforward demonstration of the result for 
Dirichlet boundary conditions mentioned above. The method, which depends on the 
result being known a priori, consists of shifting the eigenvalues by a parameter z, and 
arguing that both sides of the result to be proved are analytic in z and have the same 
zeros. Standard results from the theory of complex variables then allows the result to 
be proved. 

The first published studies devoted to obtaining formulae for functional determi- 
nants of this kind were by Levit and Smilansky [§] and Dreyfus and Dym ^U|. They 
extended the known results to more general kinds of operators, for example those of 
higher order, but continued to assume Dirichlet boundary conditions. While Dreyfus 
and Dym used a method of proof similar to Coleman's, Levit and Smilansky proved 
the result by showing that the logarithmic derivatives of both sides of the required 
result were equal. This was also the approach adopted by Forman [11] . who extended 
the known results to situations with more general boundary conditions. These results 
were rederived by yet another method a few years later [T!2j . Burghelea et al inves- 
tigated these and closely related questions using more abstract language [T3J E] (see 
also PH]); Ref JH] contains an accessible discussion of some of their results. In partic- 
ular unknown prefactors are determined in these articles. Other authors [T7]-[2D] were 
concerned with defining generalised forms of determinants (the "p-determinant" or the 
"A-determinant" ) which act as regularised versions of the conventional determinant, 
and may be useful in cases where the usual determinant does not exist. Recent work 
has extended some of this rigorous work to regular singular potentials 21 j and to more 
general kinds of boundary conditions [22] ■ Also explicit formulae that relate quotients 
of determinants with the corresponding Green function have been given 

In contrast to the many different techniques devised to prove these rather mathe- 
matically attractive results, the task of obtaining analogs for functional determinants 
with the zero mode extracted has been neglected. Apart from the paper by McKane 
and Tarlie already mentioned [3] , the only published articles to address this point have 
been by Kleinert and Chervyakov flQ [22] . They base their approach on the Wronski 
method for constructing Green functions, but they do not discuss cases where there is 
more than one variable, and again introduce an explicit regularisation as an intermedi- 
ate step in the calculation. In this article we will use the contour integration methods 
developed in [211 12Z] to avoid introducing ad hoc regularisation procedures for the zero 
modes. For a recent review of these methods see Kirsten [2H]- This will allow us to 
easily obtain the formulae for functional determinants we seek in terms of solutions 
of the associated homogeneous equation and the boundary data of the problem. The 
results generalise known answers where the absence of zero modes was assumed. 

The outline of the paper is as follows. In Section [21 we describe the method in the 
most straightforward case: obtaining the functional determinant of a second order dif- 
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ferential operator with Dirichlet boundary conditions. Having illustrated the method, 
in Section 01 we discuss how the formalism has to be modified to give results when any 
zero mode which exists has been omitted. In sections 0] and El we extend the treatment 
of sections El and El respectively, to general boundary conditions. The extension to sys- 
tems of differential equations will be briefly discussed in section El but, for simplicity, 
we will limit ourselves to an outline and give results which for the most part apply to a 
specific example. We conclude in section Two appendices contain technical results 
pertaining to the development of some aspects of the formalism described in sections 
El and El 



2 Dirichlet boundary conditions 

We begin by considering the simplest class of differential operators: 

L, = + Rjix) (1) 

on the interval / = [0, 1]. Since the determinant of any such operator will diverge, it is 
usual to calculate ratios of determinants both of which involve operators having similar 
structure, for instance, defined on the same interval and having the same leading order 
term. The index j — 1,2 labels these two determinants. We will adopt the convention 
that the j = 1 operator is the one for which we wish to calculate the determinant 
and the other, j = 2 one, will be some standard reference operator, in terms of which 
our result can be given. As long as the interval in the problem of interest is finite, 
we can shift and scale x so that x G [0, 1]. The operators in this section will have no 
first derivative term and no non-trivial function of x multiplying the second derivative 
term. We will factor out any constant terms from the Lj so that the coefficients of the 
second derivatives is —1. The functions Rj{x) will be assumed to be continuous and 
all the eigenvalues will be assumed to be positive. 

In applications which involve the evaluation of functional integrals, ratios of deter- 
minants frequently arise, as do Dirichlet boundary conditions. So it is natural that the 
first results which were obtained should have been for this case. As we have already 
mentioned in Section Q the result involves the solutions of the associated homoge- 
neous equation. In fact it is given in terms of the solution which satisfies given initial 
conditions: 

detLi = yi(l) . . 

detL 2 y 2 {\y 1 ' 

where yj(x) is the unique solution of 

L jyj (x) = 0, %(0) = 0, y' 3 (0) = l. 

The main purpose of this section is to rederive this result by a direct contour integration 
method. 
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To do this let Uj t k(x) be the unique solution of 



(Lj — k 2 )uj t k{x) = 

satisfying 

u jtk {0) = 0, u' jJe (0) = 1, (3) 

k being a complex parameter. The eigenvalues are then fixed by imposing 

= 0. (4) 

The second condition in Q is simply a choice of normalisation for the eigenfunction 
Uj,k{x); we will indicate at the end of this section what happens if we make a more 
general choice. In what follows we will assume that Wj,a;(1) is analytic in k for ^tk > 0. 
The definition of this eigenvalue problem allows the zeta function of Lj to be written 

as 

1 f ., , o„ d 



The contour 7 is counterclockwise and encloses all eigenvalues on the positive (by 
assumption) real axis. As given, the representation is valid for ^ts > 1/2. Clearly 

1 f „ , -2, d 1 u hk (l) 



-ClM = — / dMr^-ln-^-Z (5) 

The next step is to deform the contour to the imaginary axis. In the case being 
considered here there are no negative or zero modes, this means Wj,o(l) 7^ 0, and so no 
contribution arises from the origin. The leading \k\ — > 00 behaviour is independent of 
Rj(x) and is governed by 



u j>h (x) ~ sin(kx) (l + Oik- 1 )) , (6) 
j = 1, 2, so that for \k\ — >• 00 we have 

With e > a small real number introduced to indicate that the contour approaches 
the imaginary axis from the right, this allows us to write 



C Ll (s) - Cl 2 (s) = — dk(ik + e)-»^]n-*y-, (7) 

which is valid for > —1/2. 

The next relevant observation is that Uj^k{x) and Uj_ik{x) satisfy the same differ- 
ential equation, and since they are the unique solution one has Uj^x) = Uj-n-(x). 
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The integral in (JJJ) may be broken up into two contributions, k > and k < 0. The 
two contributions combine to give 

a,w-a ! w = ^r**- 2 4'n ; " u ' (1) 



7T Jo dk U 2 ,ikW 

which is valid for —1/2 < 3?s < 1/2. The upper restriction comes from the lower 
integration bound k = 0, whereas the lower restriction is from the upper integration 
limit as discussed. 

From here it is immediate that 

CUO) - CUO) = r<ik± in HMiffi = - in Sfifii = - 1„ ^, 
2 Jo dk u 2 ,ifc(l) w 2 , (l) 2/2(1) 

which reproduces eq. (J2J) using the definition 

det(X 1 L 2 - 1 ) = exp-«k(0) - Ci 2 (0)). (8) 

If we do not normalise Wj )fc (0) = 1 but instead ask that u'j k (0) = Cj, the above 
integral will give a contribution at infinity: 

CI, (0) - CUO) = In SfeW - 1„ ^mM = - ,„ c ^,o(l) 



W3,ioo(l) M 2 , (l) CiM 2 , (l)' 

because now 

^l,ioo(l) _ Ci 
M2,ioo(l) C 2 ' 

This is consistent with the answer expected: 

det(W ) = £2^11 = V|W(°) (9) 

The strength of the method is that, as discussed later, it can be applied with very little 
modification to more complicated operators and boundary conditions. However, just 
as important, we can use it to determine the analogous result to for the situation 
where a zero mode is present and has been extracted. 



3 Omitting the zero mode 

To explain how the procedure discussed in the last section can be modified to deal 
with a situation in which a zero mode is present, but has been omitted from the 
evaluation of the determinant, we stay with the simple operator (0) and Dirichlet 
boundary conditions. However, as before, we will see that the method generalises to 
more complicated situations. In addition to (J3J) we now also suppose that %,o(l) = 0, 
that is, Li has an eigenfunction with zero eigenvalue, and all the other eigenvalues 
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are assumed to be positive. The integrand in (J5J) will now have a pole at k = and 
therefore the contour 7 can no longer be deformed along the imaginary axis. Our 
approach will be to consider a modified form of the integrand which has no pole, and 
which still gives us the determinant of L\, but now with the zero eigenvalue extracted. 

We first need to determine the behaviour of iti,fc(l) for small k in order to eliminate 
the pole. To do this, we note that integrating the left-hand side of 



1 rl 

dxu\ q(x)*LiU\ kix) = k 2 / dxuio(x)*ui k(x), 

Jo 



by parts gives 



u' 10 (x)*u lik (x) - u[ tk (x)*ui }0 (x) 



1 

■ ./,) 



1 

dar«i fc (ar)(Litiio(a:))* = k 2 (u lfi \u hk ). 



Here * denotes complex conjugation and we have introduced the Hilbert space product 
( I ) on C 2 (I) by 

(u\v) — dx u(x)* v(x). 
Jo 

Although we will be mostly concerned with real functions in this paper, the general- 
isation to complex functions will be required in section El and so it is useful to allow 
for this possibility in the development of the formalism. 
Using the boundary conditions this gives 

u i,o\ l ) 

Since fi^ is finite and non-zero as k — > 0, we have the desired behaviour of t*i,fc(l). 
The minus sign has been included in the definition of f\ y k for later convenience. It is 
important to note that (JTUJl is true for any k — no small k assumption was made to 
derive it. 

We can now modify the discussion of section |2] to cover the case when a zero mode 
is present, by using the following two observations: 

(i) The function fi ik , defined by (fTU|) . vanishes at all values of k for which k 2 is a 
positive eigenvalue, as Wi,fc(l) does. Only in the case of the zero eigenvalue do we 
have the situation where fi t o 7^ 0, but = 0. Therefore if we use fi ^ instead 
of in (0), and choose the contour so as to surround the positive eigenvalues, 
the formalism developed in section |21 can be used to obtain the required result. 

(ii) The notion that Ui,fe(l) may be straightforwardly substituted by /i^ has to be 
amended slightly, since these functions have different behaviours as k —>■ 00. This 
is because the method of proof we used in the last section relies on Ui,ifc(l) and 
^2,iA:(l) having the same behaviour in this limit. This difference in behaviour is 
accounted for if we replace Ui,fc(l) not by fi^, but by (1 — k 2 )f\ t f,, since then 
Wi,ifc(l) ~ k 2 fi ik as k —>■ 00, while (1 — k 2 )fi^ remains non-zero at k — 0. 
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These remarks lead us to consider the contour integral 



. dkk- 2s ^-ln(l-k 2 )f hk , (11) 
zm J~f ak 

where the contour surrounds all of the eigenvalues on the positive k axis. This integral 
equals 

ClxW + ^t / dkk- 2s ^-\n(l-k 2 ), (12) 
2m J-y alt 

where it is understood that the zero mode has been omitted from the definition of the 
zeta-function. The second term in (|12j) is equal to 1 if the contour surrounds k — 1, and 
zero if it does not. Our final result does not depend on this choice so, for definiteness, 
we choose the contour to surround this point. 

The results (fTTj) and (|T2*|) can now be combined with the conventional result for the 
zeta function of L 2 (the eigenvalues of L 2 are assumed to be all positive). The contour 
can safely be deformed to the imaginary axis and the contribution from infinity in the 
right-half plane shown to vanish. Using fi^ = /1 and proceeding as before we find 

ClA s ) -Cl 2 {s) = / dkk —In — 1. (13) 

The conditions on s for this representation to be valid are the same as those given in 
section |21 and therefore we may differentiate with respect to s and set s = 0. If we 
impose the condition Mj jfc (0) = 1, then the argument of the logarithm in (fT3*j) equals 
Wi,ioo(l)/w2,«oo(l) = 1 in the limit k — > 00. Therefore the contribution from the upper 
limit of integration is zero. So only the contribution from the lower limit remains and 
we obtain 

CU0)-CU0) = -ln-^. 

Using (JBJ and (fTUjl this gives 

det'Li / 10 <2/i |2/i> 



detL 2 w 2i0 (l) yi(l)*y 2 (l) 



(14) 



where we have denoted the determinant of L\ with the zero mode extracted by det'la. 

This result should be compared to (J2J), rather than Q, since we have imposed the 
normalisation condition y'j(0) = 1. If we had not imposed normalisation conditions, 
an extra factor of Ci/c 2 would be present as in (jUJ). We see that the modification 
required when a zero mode is extracted gives a final result which is still very simple, 
involving only the derivative of the zero mode at an end-point (both end-points if the 
function is not normalised) and its norm. It agrees with the result derived in Ref. [5], 
once typographical errors in eqns. (2.7) and (2.14) of that paper have been corrected. 
However, when making the comparison it should be borne in mind that the operators 
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used in that paper are minus those used here, and so the extraction of a mode gives a 
result which differs from ours by a sign. We can check that our result has the correct 
sign by noting that the lowest eigenfunction should have no nodes, so y[(0) and y[(l) 
should have opposite signs. In particular, if we impose the condition y[(0) = 1, then 
1/1(1) < 0. Therefore, (fTljl should be positive, which is correct, since all eigenvalues are 
positive. 

In applications it may not even be necessary to calculate (yi\yi). This is because, 
if the zero mode has been extracted using the method of collective coordinates [3], a 
Jacobian appears due to the transformation to new coordinates which, to lowest order, 
is proportional to the norm of the zero mode, and cancels it out. A little care is required 
with this cancellation, since the zero mode appearing in the Jacobian function will have 
a definite normalisation associated with it, which it inherited from the classical solution 
(for instance, the zero mode will be the spatial derivative of the classical solution if 
translational invariance has been broken). Therefore an overall constant may have to 
be factored out of the zero mode to get yi, and similarly from the Jacobian to get 
(Vi\yi), but these are purely algebraic procedures — no integration will be required. 

To illustrate the ideas presented so far in this paper, let us give an example. A 
rich source of operators of the type (JIJ are Hamiltonian operators in one-dimensional 
quantum mechanics. The functions Rj{x) will be the potentials; we will take as our 
example the linear potential jHOj Ri{x) = x — xq and normalise it with respect to the 
case R2{x) = 0. This potential has been chosen purely on grounds of simplicity; it 
has not to our knowledge occurred in a calculation of fluctuations about a nontrivial 
solution as described in the Introduction. As a consequence, if we wish the lowest 
eigenvalue to be zero, the constant x has to be appropriately tuned. 

The general solution of the equation u'[ = (x — x — k 2 )u\ which satisfies «i,fc(0) = 
is, up to normalisation, 

ui,k{x) = Ai(—z) Bi(x — z) — Bi(— z) Ai(x — z), z = x + k 2 , (15) 

where Ai and Bi are Airy functions [3*T] . Implementation of the second boundary 
condition, «x,fc(l) = 0, fixes the condition for k 2 to be an eigenvalue. This condition 
takes the form of a transcendental equation involving Airy functions, however a good 
approximation to the allowed values of k 2 may be obtained by using the asymptotic 
form for Airy functions, valid when the argument of the function is large |3lJ. One 
finds z 3 / 2 = (z — I) 3 / 2 + 3mr/2, n = 1,2, . . ., which implies that z « n 2 7r 2 + 1/2. 
Although this result has been formally derived under the assumption that z is large, 
it turns out to be an excellent approximation even for quite small n: correct to 0.01%. 
For the n = 1 case, which is the least accurate, the true result is 10.3685 . . . and the 
above approximation gives it 2 + 0.5 = 10.3696 . . .. Note that for large n the eigenvalue 
condition becomes z ~ n 2 TT 2 , which is the result which would be found in the case 
R2{x) = 0. Thus it is at least plausible that the ratio of the product of eigenvalues 
converges. However, the whole point of the formalism we have been developing is to 
bypass the determination of the eigenvalues, and simply work with the solutions of the 
corresponding homogeneous equation. 
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The relevant solution of the homogeneous equation for this particular example is 
Ui{x) = lim fc ^ Ui t k(x), where Wi,fc(x) is given by (fT5j). but subject to the normalisation 
condition u' lk (0) = 1, as in (JHJ). Thus 

Ai(— x ) Bi(x — Xq) — Bi(— x ) Ai(x — x ) 

Viix) = 



Ai(-xo) Bi'(-x ) - Bi(-x ) Ai'(-x ) 
7i [Ai(— xo)Bi(x — xq) — Bi(— xo)Ai(x — xq)] , (16) 



where we used the result for the Wronskian determinant [31 



W {Ai(z),Bi(z)} = -. 

71 



Similarly, 

Therefore using (jSJ), 

det Li 
det L 2 

If xo = 0, the result simplifies considerably |TTTj : 
detLi T(l/3) 



y 2 (x)=x. (17) 



7T [Ai(-x )Bi(l - x ) - Bi(-x )Ai(l - x )] . (18) 



det L 2 2 3V6 



Bi(l) - v / 3Ai(l) = 1.085. (19) 



Now suppose that we artificially fix xo to be Xq = 10.3685 ... — the smallest value 
of z for which u\^{x) given by (fT5j) satisfies Wi,fc(l) = 0. Then k 2 = is an eigenvalue, 
yi(l) = and the expression in ()18|) is identically zero. To find the ratio of functional 
determinants with the zero mode extracted from the determinant in the numerator, we 
use (P3J) and (IEJ) to find 

¥^ = -^T^^^ ¥^7 VTV u = 0-050666. (20) 

detL 2 7r(Ai(-xS)Bi'(l - arg) - Bi(-x^)Ai'(l - x* )) v ; 

Let us mention that (yi\yx) can be expressed as a combination of Airy functions and 
their derivatives. However, no simplifications occur and we do not display the resulting 
answer. 



4 General boundary conditions 

In the last two sections we have illustrated how an expression for the ratio of functional 
determinants can be derived in the simple case of a single differential equation of the 
type (0) and Dirichlet boundary conditions. In the remainder of the paper we will 
extend this treatment to more general situations. Those parts of the proof involving 
zeta-f unctions and contour integration will remain essentially unchanged; the novel 
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aspects arise from setting up the problem so that the method naturally generalises. 
Since the main aim of this paper is to provide a transparent discussion available to 
a wide readership, we will limit our discussion to more general boundary discussions, 
leaving the topic of more complicated second order operators to a future paper We 
will, however, briefly discuss a system of two coupled operators of type (JTJ) in section® 
in order to provide some indication of how the procedure generalises in that direction. 

The boundary conditions on the functions on which Lj operates fall into two classes: 
separated and non-separated. Separated boundary conditions are defined as those 
which do not mix conditions on Uj tk (x) and u'j k (x) at the boundary x = with con- 
ditions at the boundary x = 1. Robin boundary conditions Auj tk (0) + Bu f j k (0) = 
0, Cti Jt fc(l) + Du'j k (l) = 0, where A,B,C,D are given real constants, are the generic 
case. Non-separated conditions can involve Uj jk (x) and its derivative at different bound- 
aries in the same equation. Periodic boundary conditions Uj jk (0) = Uj ik (l),Uj k (0) = 
u'j k {l) are a common example. In previous sections, where Dirichlet boundary con- 
ditions were imposed, initial conditions were specified at x = which guaranteed a 
unique Uj tk (x) on which to impose boundary conditions. The imposition of more com- 
plicated types of boundary conditions, and in particular those of the non-separated 
type, does not permit us to proceed in the same way, and we need to set up a more 
systematic approach. 

The first modification consists of no longer working with the second order differential 
equation 

d 2 

= (Rj(x) - k 2 ) Uj , k , (21) 

but going over to the first order formalism where we define Vj )k (x) = u'j k (x), and view 
the column vector 

■**>-(:$)) (22) 

as the basic element of the theory. Then (}2T)) may be written as 



d_fu jtk (x)\ _( 1 \ fu j:k (x)\ 
dx {v j>k (x)J ~\R 3 -k 2 ) {v j>k (x)J- 



(23) 



The general form for the boundary conditions can then be expressed in terms of Uj tk (x), 
using the notation of Ref . [TT] , as 

«M)+»M1 = (1. (24) 

where M and iV are 2x2 matrices. For example, in the case of Dirichlet boundary 
conditions 
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Now let ufl(x) and ufl(x) be two independent solutions of the differential equation 
We will not impose the boundary condition (|2"^) on these solutions, instead they 
will be uniquely specified by the imposition of initial conditions. They may then be 
thought of as basis functions from which a general solution can be constructed: 

f Uj , k (x)\ = (aufl(x) + /3v®{x)\ 
\v j>k ( x )) {av%(x)+Pv$(x)) 

Wf \ ( 2 ), 




(26) 



In the special case when k 2 is an eigenvalue, the constants a and (3 can then be chosen 
so that Uj jk (x) satisfies the boundary conditions (J2IJ). 
Now 

= | it] "111 I (27) 
is a fundamental solution to (|2^|) . Since det Hj jk (x) ^ 



So, in summary, 



^(X)J-^ (X) ^ (0) U(0), 



is the unique solution of (|2^j) with initial value (uj tk (0),Vj tk (0)) at x = 0. This allows 
us to rewrite the boundary condition (J23|) as 

(M + NH hk {l)Hrlm M°j) = 0. (28) 

Nontrivial solutions only exist if 

det (M + NH j>k (l)H-£{0)) = 0, (29) 

and this gives the condition for k 2 to be an eigenvalue. 

Although we could in principle work with any set of {u^ k \x), uf k \x)}, a particularly 
suitable choice is 

$8 )-(-)■ 

This has a number of simplifying features. For instance, det Hj >k (x) = 1 for all x, 
a — M i,fc(0),/3 = ^,A;(0) and the eigenvalue condition becomes det(M + NHj ik (l)) = 0. 
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A simple example of a set of functions satisfying (j3~Uj) are those for the equation with 
Rj = 0. They are given by u^ k (x) = coskx and u^l(x) = k~ x sin kx. 

The solution Uj fc (x), by contrast, depends on the nature of the boundary conditions. 
For Robin boundary conditions, for example, Aa + Bf3 = 0, and we may take a = —B 
and (3 = A, Uj ;k (x) being only defined up to a constant as usual. Thus 

u jk (x) = -Bu ( *l(x)+Aufl(x). (31) 

It is now straightforward to calculate det(M+NHj tk (l)) using the matrices appropriate 
for Robin conditions: 

M ={ A c n). »=(° n). (32) 



V / 

to find that 

det(M + NH jik (l)) = Cu jtk (l) + Dv jjk (l), 

so that imposing the boundary conditions at the "final" point gives the eigenvalue 
condition (|25jl. 

For non-separated boundary conditions, initial and final conditions are mixed to- 
gether and the appropriate linear combination (j2Bj) is determined by imposing either 
one of the boundary conditions. For example in the case of periodic boundary condi- 
tions, Uj jk (0) = Uj jk (l) implies that at(l — «^(1)) = f3u^ k (l), so that we may take 

ujM = ^S«4^) + (i - *43(i)) u< il( x )- ( 33 ) 

The other boundary condition then gives the eigenvalue condition. Explicitly calcu- 
lating det(M + NH jtk (l)) using O and M = -N = J 2 , where J 2 is the 2 x 2 unit 
matrix, one finds 

det(M + NH jik (l)) = v jjk (0) - v jik (l). 

It is now possible to see how the approach adopted to investigate the case of Dirichlet 
boundary conditions fits into this more general structure. The solution Uj^{x) which 
satisfies (J2l?j) is simply equal to u^l(x) introduced in this section, that is, a — 0, /3 — 1. 

Using (|23j) . det(M + NHj !k (l)) = it) 1(1), so that the eigenvalues are fixed by imposing 
Uj t k(l) = 0, just as in (J3J). For Neumann boundary conditions a = 1,(3 = 0; for other 
cases both of the fundamental solutions u^l(x) and u^l(x) will be required to construct 
a solution Uj tk (x) which has the correct properties. 

The discussion presented in this section so far indicates how our earlier proof can 
be extended to more general types of boundary conditions. Assuming suitable analyt- 
icity properties in k, the eigenvalue condition (0J) is now replaced by the more general 
condition (j2*U|) . and so the analogue of (0) becomes 

> /n > /x 1 f „ r 2 , d i det(Af + JVfT lfc (l)) 
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We next perform the same steps as in section 2. The leading \k\ — > oo behaviour of 
u^l(x) is governed by eq. (0), whereas for we have 

ufl{x) ~ cos(fcx)(l + 0(k- 1 )). 

This shows 

d_ det(M + NH lik (l)) 2 
dfc det(M + jV# 2ifc (l)) 1 

and all goes through as before. Thus, with the definition 

lim Ha k (x) 
yf [x) ) fc-o ' 

the formalism immediately provides the answer for the quotient of determinants 

detLi _ det(M + iVr 1 (l)) 
det L 2 ~ det(M + NY 2 (1))' 

From a practical point of view, in order to calculate (|35|) only the matrices M and 
N, which specify the boundary conditions, and the two independent solutions {yf\x), 
yf\x)} of the homogeneous equation Ljyj = are needed. These fundamental solu- 
tions also satisfy the initial conditions (J3Uj) inherited from the solutions {u^l(x), u °l(x)}: 

m - { y f"(o) J " { o i j ■ (36) 

For instance, in the example introduced in section 3 where Ri (x) = x, the fundamental 
solutions are 

yPfr) = [Bi(ar) + v/3Ai(x) 



(35) 



2/i 0*0 



£(1/3) 
2 3V6 



Bi(z) - v^Ai^)! . (37) 



For Dirichlet boundary conditions y±(x) is simply y^\x) (compare with (jl9|0 . and for 
Neumann boundary conditions it is y[ (x). 

Let us end this section by giving an alternative, slightly simpler, form for 
which will be useful in the next section. We have seen in the case of Robin and 
periodic boundary conditions, that if we implement one of the boundary conditions, 
then det(M+iViifj j fc(l)) is proportional (equal, with the correct choice of normalisation) 
to the other condition. This will be true in general. To explicitly show this, let us write 
the condition (EH) out in full: 



m fl2U) out in full: 

m n u j; k{0) +m 12 v j:k (0) +n u u j:k (l) + n 12 v j>k (l) \ _ / \ . 
m 21 u j:k (0) + m 22 v jik (0) + n 2X u^ k (l) + n 22 v jik (l) J V / 
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We ask that satisfy only one boundary condition. This could be either one, but 

one choice gives 

mnu jtk (0) + mi 2 fj,fc(0) + n u u jtk (l) + n 12 v jik (l) = 0. (39) 

We now want to construct ) which explicitly satisfies the condition (}3*9*j) . To do 

so we expand Uj tk (x) as in (J26|) and implement (J39j) . This gives (with an appropriate 
choice of normalisation) 



mi2 + n u u ( ll(l) + n 12 vf k (1)] uf}(x) 
mil + nu«$(l) + ^SW] «$(*)■ ( 4 °) 



+ 



Using the expression ()4()j) to calculate an explicit expression for the second row of the 
left-hand side of (JHBJ) gives 

detM + detiV + iiv](l)(num 2 2 — n 2 iWi 2 ) + ttp(l)(mun 2 i — m 2 in 11 ) 
+ vfj,(l){n 12 m 22 - n 22 m 12 ) + v^ ] k (l)(m u n 22 - m 21 n 12 ), 

which is simply det(M + NHj yk (l)) expanded out and simplified using det Hj tk (l) = 1- 
Therefore we have shown that 

det(M + NH jtk (l)) = m 21 u j)k (®) + m 22 v jyk {0) + n 21 u jik {l) + n 22 v jik {l). (41) 

The fact that the constant of proportionality between the two sides of this equation is 
1 is a consequence of the judicious choice of normalisation of f!40|) . If we had chosen 
to implement the second row of (}3*Hj) instead of ()39|) . we would simply have obtained a 
different Uj tk (x), which, when substituted into the first row of (JHHJ), would have again 
given det(M + NH j<k {\)). 

As a result of (j4~Tj) we find for the quotient of functional determinants 

detLi _ m 2 ig/i(0) + m 22 ^(0) + ngiyiQj + n 22^(l) 
detL 2 m 21 y 2 (0) + m 22 y' 2 (0) +n 21 y 2 (l) + n 22 y 2 (l)' 

which is the alternative form to ()35)1 mentioned above. It is now possible to immediately 
write down the required formulae for any boundary condition. For instance, in the case 
of Robin boundary conditions defined by (JH2I), 

detL x = C yi (l) + Dy[(l) 

detL 2 Cy 2 (l)+Dy' 2 {iy 1 ' 

and for periodic boundary conditions defined by M = —N = I 2 , 

detLx y'M-y'^l) 



detL 2 ^(0)-^(l)" 



(44) 
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5 Omitting the zero mode: general case 



In this section we will derive the result for the ratio of functional determinants when a 
zero mode has been extracted, in the case when general boundary conditions of the type 
(jUj) have been imposed on operators of the type (jTJ). As before, we will assume that 
the operator Lj with the boundary conditions imposed defines a self-adjoint operator. 
This requires that the elements of the matrices M and iV satisfy certain conditions, as 
described in Appendix A. The main result which we will require will be the nature of the 
generalisation of (fTU|) . where ui,fc(l) will presumably be replaced by det(M + NHj tk (l)). 
As in section |3J the starting point of the analysis is 

[v ly0 (x)*ui, k (x) - v lt k(x)ui i0 (xy]l = k 2 («i, |wi,fc) • (45) 

We implement both boundary conditions on the k = solution, but only one of the 
boundary conditions on the solution for general k, since we wish to keep k arbitrary, 
and not restrict it so that k 2 is an eigenvalue. 

To gain some insight, we first look at particular boundary conditions before we 
outline the procedure for general boundary conditions. For Robin boundary conditions 
we have Au lt0 (0) + Bv lfi (Q) = 0, Cu lfl (l) + Dv lt0 (l) = and Au hk (0) + Bv hk (0) = 0. 
The left-hand side of f!45|) then gives zero for the contribution at x = (both functions 
satisfy the boundary conditions there) and at x — 1 gives 

- ^ 1)0 (i)x*(i) - ^(iKo(i)* = - (^r-) [cVfcW + Dv iM] ■ 

Therefore 

Cu 1>k {\) + Dv 1)k {l) = - ( ° f .A k 2 ( Ul , \u hk ). (46) 

V M i,ol-U / 

For periodic boundary conditions we have ui,o(0) = ^1^(1), fi,o(0) = t>i 5 o(l) and 
^i,fc(0) = iti,fc(l). The left-hand side of (|43j) then gives 

- ui,fc(lHo(0)* + ui,fc(0)u ll0 (0)* 

which leads to 

vi >k (0) - v 1>k (l) = ^— fc 2 (« 1|Q |'Ui ifc ). (47) 

For all special examples considered, the result takes the form (see eqs. ([TO]). (jUj)) 

det(M + NH hk (l)) = Bk 2 (u lfi M } (48) 

where B depends on the zero mode data at the boundary and on the boundary condition 
imposed. 

In fact, the procedure described above can be extended to the case of general 
boundary conditions and the analogous results will always have the form (J48|) . To see 
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this recall from the last section that implementing one boundary condition only leads, 
in the case of general boundary conditions, to (see (J3~9*)) and (jHJ)) 



mu«i,fc(0) +m 12 vi )fc (0) +nn«i )fc (l) +ni 2 vi,fe(l) = (49) 
m 2 i«i,fc(0) +m 2 2«i,*(0)+n 2 i« ll fc(l) +n22«i,fc(l) = det(M + NH 1>k (l)). 

When k — 0, Ui,fc(x) is the zero mode and so (|4*9^1 holds, but with det(M + NH 10 (1)) 
replaced by 0. 

We wish to obtain expressions for ui tk (x), ui t o(x), and their derivatives, at the 
boundaries in order to substitute them into the left-hand side of and ultimately 
obtain an expression for det(M + NH l k (l)). To do this we have to proceed by cases, 
which depend on the particular boundary conditions prescribed. Suppose that, for 
instance, the boundary conditions are such that vn\ 2 n 22 — ni 2 m 22 7^ 0. Then in this 
case we rewrite (|4*H|) as 




«i,k(0)\ = / \_(m n n n \ (uiA^Y 

v hk (l)J {det(M + NH 1:k (l))J { m 21 n 21 J {u 1 , k (l) / 



(50) 



The condition on the elements of M and N that we have imposed implies that the 
matrix on the extreme left of (J50|) is invertible and so we may solve for Vi ik (0) and 
fi 5 fc(l). These intermediate results are of no intrinsic interest and are given in the 
Appendix A. Obviously the equations for fi,o(0) and t>i,o(l) are identical to those for 
general k, except that det(M + NH 1>k (l)) is replaced by det(M + NH lfl (l)) = 0. We 
may now eliminate Vi ;k (0) , v 1)k (l) , Vi >0 (0) and t>i,o(l) from using these equations. 
The basic steps are outlined in Appendix A. We obtain 

det(M + NH l k (l)) = ^12^22 m 12 n 2 2 ( M Q | M ^ n i2 m 22 - m i2 n 22 ^ 0.(51) 

mi 2 «i, (l)* + n 12 u lfi {0)* 

Expressing {v i,fc(0), v i,fc(l)} in terms of {mi,a;(0), ^1^(1)} is just one possibility; there are 
five more ways to express two boundary data by the complementary ones. Proceeding 
as described gives the list of results presented in Appendix A. It is seen that the precise 
form that B takes depends on which combination of elements of M and iV are non-zero. 
The formulae in Appendix A, together with (|51|). cover all possible cases. If more than 
one of the results is appropriate, then they are, of course, equivalent. 

Having established that (}4*5j) holds in general, we may simply replace (fTU|) by 

det(M + NH 1>k (l)) = -k 2 f hk ; f 1>k = -B{u lfi \u lik ) , (52) 

and repeat the steps in sectional One finds 



CU0)-CU0) = -ln 



/i,o 



det(M + iVF 2 (l))' 
which using (JHJ gives 

det' Li B (y\\y\) 

detL 2 " det(M + NY 2 (1)) ' 



(53) 
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where B takes on the following values: 



n 12 m 22 - m 12 n 22 

771122/1 (1)* + 71122/1 (0)*' 

m n n 2 i - m 21 n u 

mu2/i(l)* + nu2/i(0)*' 
m X \vn 22 - mi 2 m 21 

771122/1(0)* + mn2/i(0)*' 
m 12 n 21 - m 22 nii 

777122/U 1 )* - 77112/1 (0)*' 
7*12^21 - 7l 2 277ll 

77 1 22/U 1 )* + 77ll2/l(l)* 

777nn 22 - m 21 n 12 

77122/U )* - 777,ii2/i(l)* ' 



it n 12 m 22 


~ 777l 2 7l22 




0. 


if mnn 2 i 


- 771217711 




0. 


if mnm 2 2 - 


- m 12 m 21 




0. 


if m 12 n 2 i 


- m 22 n u 




0. 


if ni 2 n 2 i 


- n 22 n u 




0. 


if m u n 22 


- m 21 n 12 




0. 



(54) 



As an example, for periodic boundary conditions when M = I 2 , the third of the 
above results is applicable, and therefore B = 1/2/1(0)*. Equivalently the fifth result 
may be used, since N = — J 2 , to give B = 1/2/1(1)*. Furthermore det(M + NY 2 (1)) = 

2/2( ) - 2/2(1) ( see (El); so that 

det' Li (2/1 1 2/i) f ,,v 

detL 2 2/i(0)* [2/2(0) -2/2(1)]' 1 ] 

Let us end this section by showing that (|55|) is in agreement with the result obtained 
by a different method 0. In order to show that the two results are identical, we need 
to establish a correspondence between the solutions utilised in the two approaches. 
In this paper we have taken u^l(x) and u?l(x) — which in the k — ► limit become 
2/ 1 ) (x) and y^ 2 '{x) — as the two independent solutions (to facilitate comparisons with 
the formulae given in Ref. [Hj, we will drop the j subscript for the rest of this section). 
In Ref. |H] the two independent solutions were denoted by 2/1 (x) and y 2 (x), but we will 
denote them here by y^(x) and y^ n \x). The solution y®(x) is simply what we have 
called y(x) in this paper, that is, the solution which is the zero mode if one exists, 
and which satisfies only one of the boundary conditions, if a zero mode does not exist. 
Thus from (jSHJ), when k — > 0, 



y®(x) yW(x) \_( y^{x) j/ (2) 
yV'tx) yW(x) j ~ [ y^'(x) 




(56) 



where y^(x) = jy^(x) + Sy^(x) is any solution which is independent of y^(x). 
Since detY(x) = 1 in our formalism, the Wronskian of the two solutions y®(x) and 
y^ ll \x) is 



(I) (x)2/ (II) '(x) - y {ll \x)y {l) \x) = aS - /3 7 . (57) 



All of this holds whatever the nature of the boundary conditions — they will corre- 
spond to particular choices for a and f3. Specifically, for periodic boundary conditions, 
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shows that a = y^(l) and (3 = 1— y^'(l), so that the Wronskian, W, is given 
by W = Sy {2) (l) - 7 (l - lt is straightforward to show that - y {u) {0) 

is also equal to this quantity. So, for periodic boundary conditions, and y®(x) given 
by the k -> version of (J33J), W = y (II) (l) - y (II) (0). This verifies that eq. (5.2) of 
Ref. [Sj reduces to ({HE} above, up to a minus sign which originates from an overall sign 
difference in the definition of the operators Lj. 



6 Systems of differential operators 

As we have stressed several times already, our main aim in this paper is to present 
results and proofs in such a way as to be as free of technical details as possible. So, 
while the results for a second order differential operator of the type (0) presented in 
the earlier sections of the paper generalise to systems of differential operators, we will 
only give a brief outline of the formalism here, leaving the details to a forthcoming 
publication [3]. In particular, we will, for definiteness, discuss most of the results in 
the context of a specific example, encountered in the study of transition rates between 
metastable states in superconducting rings [321 EH]- The differential operator in this 
problem is defined on the interval [—1/2, 1/2] and has the form 

;i - ^ 2 )e 2 ^ \ d 2 

dx 2 




L 1= \ \=-I 2 — + Q 1 (x). (5f 



dx 2 



+ (1 - 2tf 



In the actual problem this example is taken from, // is a wavevector, which is small 
(so that fi 2 < 1), and I is the circumference of the ring divided by the temperature- 
dependent correlation length. The solutions of (Lj — k 2 ) Uj t k(x) = will be denoted 
explicitly by 

^ w =&£!)■ (59) 

The sole purpose of the operator L 2 is to render the ratio detLi/ det L 2 finite, and it 
is chosen on grounds of simplicity and convenience. 

In the completely general case u J)fc (x) will have r components Uj t k, a (x); a = 1, . . . , r, 
and Lj will take the form of an r x r matrix containing second order derivatives along 
the diagonal. It is again convenient to go over to the equivalent problem containing 
only first order derivatives, which will have the form (|23jl . but with the column vectors 
having entries Uj t k,a{x) and Vj } k, a {x); a = 1, . . . ,r, and with the square matrix having 
the structure 

/ 



v Qj(x) - k 2 I r 



(60) 
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where I r is the r x r unit matrix and Qj{x) contains the non- derivative terms in Lj, 



as in our specific example (J58J). The boundary conditions can again be expressed as 
(I24J) . but now with M and N being 2r x 2r matrices. The advantage of defining the 
structure in this way is that much of the formalism described in section 0] goes through 
unchanged. For instance, the fundamental solution consists of 2r independent solutions 
of Li and their derivatives: 



H jM (x) 



■til 


[x) ... 


••• "jS^ 


J j,k,r 


[x) ... 


• • • U j,k,ri X 




» • • • 


■■■ V j£l( X 




» ... 


■■■ V 3\£r( X 



(61) 



and the eigenvalue condition is again (J29|) . We again choose the solutions so that 
Hj,k(—l/2) = h, just as we did in the one-component case (|3U[). When no zero modes 
are present, the representation (}3"4*j) then leads to the final result ()35|) as before. To cal- 
culate the right-hand side of we only need solutions of the homogeneous equations 
Ljjj = and the matrices M and N. 

In the particular example we are considering here, we will impose "twisted boundary 
conditions" (321 ESJ M = -diag {e^, e~ ifd , e"^} and = J 4 . It is easy to check 
that 




(62) 



is an eigenfunction of the problem with zero eigenvalue and so det(M+iV.£/i i o(//2)) = 0. 

Since a zero mode is present we have to proceed in the general way indicated in 
section but suitably generalised to this particular two-component case. The analog 
of (SSI) is 



e ko, a < 



X)U 1Xa {x) - V likta [x)U 10 AX) 



a=l 



1/2 
-1/2 



fc 2 (Ui >0 |Ui >fc ). 



(63) 



Now we assume that u 10 (a;) satisfies all of the boundary conditions, so it is a zero 
mode, but ui^{x) satisfies only 3 of the 4 — for instance, those used in Appendix B, 
see (IB.IJ) . This means that k is not restricted to values for which k 2 is an eigenvalue. 
Substituting in these conditions, (JfiHj) becomes 



«l,fc,2(«/2) 



e-^ l v 1A2 {-l/2) u* 1>0>2 {1/2) = fc 2 (u li0 |u 1)fc ). 



(64) 



Only the component of Ui^x) that does not satisfy the boundary condition survives 
on the left-hand side of the equation. If u x ^(x) is appropriately normalised, the term 
in square brackets on the left-hand side of (J64j) is det(M + NHi^(l/2)), and we find 

k 2 



det(M + NH ljh (l/2)) 



1,0,2 



(uio|ui fc ) = -k f hk 



(65) 
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where fi^ = — ^(u^olui^) and B = — | y$ i2 (l/2) } • The only remaining task is the 
calculation of the appropriately normalised function yi(x). The method used to obtain 
it is outlined in Appendix B. Using this procedure we find 



4/(1 -/i 2 )sinh 2 e il ^ 2 



»w= ^ — {. e -4 («) 

which should be contrasted with the unnormalised result (|62j) . Therefore 



(yi|yi) = 8/ 2 (l-/i 2 ) sinh 2 (|) 

yunm 



/1.0 = -TTTTTTT = 772 — > ( 67 ) 



which agrees with the result calculated in Ref. [32] using a different method. This 
result gives det'Li, provided that we divide by detL2. The most convenient choice for 
L2 is to take it to be of the form L2 = —^(cP/dx 2 ). The calculation of det L2 then 
decouples into the calculation of the determinant of two identical operators of the kind 
studied in the earlier sections of this paper. 



7 Conclusion 

The two main objectives of this paper have been firstly to present a simple, yet rig- 
orous, way of determining the now well-established formulae for ratios of functional 
determinants, and secondly to extend the results to the situation where a zero mode 
exists, but is omitted, in one of the determinants. We have stressed clarity, preferring 
to study the simplest types of operators, and leaving the most general case to a later 
publication jlj. 

The main idea behind the method of proof which we have adopted involved defining 
two solutions of the eigenvalue equation LjUj^{x) = k 2 Uj^(x), denoted by u^l(x) and 

(2) 

Ujl(x). Here L\ is the second order differential operator of interest, and L2 is a similar, 

typically simpler, operator used for normalisation. Although u^l (x) and ufj. (x) satisfy 
the eigenvalue equation, they are not, in general, eigenfunctions since they do not satisfy 
the boundary conditions (simple boundary conditions such as Dirichlet and Neumann 
are exceptions). Instead we impose "initial" conditions on these two solutions - 
conditions on the functions and their first derivatives at x = — so that they are 
uniquely defined. The particular choice we made in this paper was (|3U|). but obviously 
other choices are permissible. 

Once two unique and independent solutions of the eigenvalue equation have been 
defined, any other solution may be expressed as a linear combination of these by spec- 
ifying the two constants a and (3, as in (J26)) . In particular, if we wish to determine 
eigenfunctions we would first apply one of the boundary conditions to find the ratio 
of a to (3. Since the overall normalisation of eigenfunctions is irrelevant, this appli- 
cation of only one of the boundary conditions determines the functional form of the 
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eigenfunction Uj t k{x). The imposition of the second boundary condition has the effect 
of restricting values of k 2 to those which are eigenvalues — which are real, since the 
operators considered here are all taken to be self-adjoint. 

The proof of the results begins from the observation that the quantity det(M + 
NHj t k(l)) is zero only when k 2 is an eigenvalue. It follows that its logarithmic derivative 
has poles with unit residue at values of k such that k 2 is an eigenvalue. Therefore a 
zeta-function can be defined through a contour integration, where the contour is chosen 
to surround these values of k. Since these poles are restricted to the non-negative real 
axis, the contour may be deformed to lie on the imaginary k axis, and then evaluated 
to give the required result (|35j) . The only other information which is required during 
this process relates to the nature of the solutions of the differential equation for large 
\k\ 111. 

An alternative result, which is slightly more explicit than ()35|) . can be obtained 
by using the solution t/j(x) = lim fc ^ u j,k(x) rather than yf\x) = \im k _, u^l(x) and 
yf\x) = lim fe ^ (x) , to express the result. Here one has to be careful with questions 
of normalisation. The solution Uj yk {x) satisfies only one of the boundary conditions, but 
is unnormalised. However if its normalisation is chosen appropriately, then det(M + 
NHj^iX)) can be made to exactly equal the second boundary condition (that is, with 
a constant of proportionality equal to 1). The effect of this, after the k — > limit 
has been taken, is to allow det(M + NHj t0 (l)) in (jHSJ) to be replaced by the second 
boundary condition, where the function appearing in this boundary condition is simply 
Uj(x), appropriately normalised. In the case of the generalised boundary conditions 
considered in section HJ the appropriately normalised Uj y k{x) satisfying the first part 
of the boundary conditions is given by (J4(Jj) . the expression of det(M + NHj t k(l)) in 
terms of the second boundary condition is given by (|41|). and the final result, in which 
the ratios of functional determinants are expressed in terms of Vj(x), is (|4^|) . 

One of the advantages of the above method of proof is that it can be relatively simply 
adapted to the case when k = is an eigenvalue of L\ which is excluded from the deter- 
minant. This is achieved by noting that in these cases det(M + NHi k{l)) = —k 2 f± } k, 
where fia 7^ 0. Since /i^ will, like det(M + NHi^(l)), vanish at all the non-zero 
eigenvalues, the proof we have described can be repeated but with det(M + NHi k (l)) 
replaced by fi^ — or rather, by (1 — k 2 )fi tk , so that the large \k\ behaviour remains 
unchanged. The result will be a formula for the ratio of functional determinants, but 
with the zero mode omitted in the determinant in the numerator. It is given by (J53|) . 
and although it depends on a constant B whose form depends on the nature of the 
boundary conditions, it is nonetheless simple, in the sense that it only depends on (i) 
the zero mode yi(x), and (ii) the matrices M and N which define the boundary condi- 
tions. Furthermore, we expect that in many applications the norm (yi\yi) will cancel 
with a similar factor resulting from the transformation to a set of coordinates which 
allows easy extraction of the zero mode. So, once again, only the values of solutions of 
the homogeneous equation at the boundary are required. 

In the examples given in this paper, we have been able to obtain the solutions of 
the homogeneous equations analytically, and so determine their values, and/or those of 
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their derivatives, at the boundaries. In most cases of interest this will not be possible, 
but this does not create any real difficulty: it is straightforward enough to determine 
the solutions yf\x) and yf\x) numerically and so calculate the required result from 
either or Since functional determinants of this sort are often perceived 

as being the most difficult quantities to evaluate in calculations of fluctuations about 
classical field configurations, we hope that the results presented in this paper will prove 
useful in this context. 
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A Determinants in the presence of zero modes for 
general boundary conditions 

In this section we give details of some of the algebraic manipulations used in section 
and collect together the results for the different cases. 

Let us first consider the case discussed in the main part of the text, that is, boundary 
conditions which are such that m\ 2 n 22 —m 22 ri\ 2 ^ 0. This condition allows us to rewrite 
P2D as 

Vfc(0)\ ( m 12 n 12 Y 1 ( \ 

^ V m 22 n 22 ) {det(M + NH 1>k (l))) 

m 21 n 21 J \u 1>k (l))' 
Multiplying this out yields 

v hk (0) = f {-n 12 det(M + NH ltk (l)) 

[m 12 n 22 - n 12 m 22 \ 

- [mun 22 - n 12 m 21 ) u 1>k (0) - [nun 22 - n 12 n 21 ] u 1;k (l)} 

and 

vi, k {l) = 7 r {mi2 det(M + NH lk (i)) 

[m 12 n 22 - n l2 m 22 \ 

- [m 12 m 21 - m u m 22 ] u lik (0) - [m 12 n 21 - n n %] «i,fc(l)} • 

Two similar equations hold for t>i,o(0) and Vi j0 (l), but with det(M + N H l k (l)) replaced 
by 0. Using these four equations we find 

r r \* / \ r \ ( [^i2^i,o(l)* + wi 2 Mi,o(0)*] , n ATU , 

[vifl(x) u lik (x) - v ljk (x)u lfi {x) L = = = det(M + NH ljk (l)) 

[m 12 n 22 - n l2 m 22 \ 

detiV (detM) N 



+Ml,fc(l)«l,o(0)* 

+ui,*(0)ui,o(l)'' 

+u 1>k (o)u lfi (oy 

If k is an eigenvalue, the first term in the above expression vanishes. Therefore, in the 
case that k is an eigenvalue, in order to have a self-adjoint problem the sum of the last 
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m 12 Tl 22 
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(m 12 n 22 - n 12 m 22 )* 


m 12 n 21 - 


num 22 


(m 12 n 21 -nnTO 22 )*| 


m 12 n 22 - 


ni 2 m 22 


(m 12 n 22 -n 12 m 22 )* J 


(mun 22 - 


- n 12 m 21 )* 


(m u n 22 - n 12 m 2 i) 


{m l2 n 22 - 


- n 12 m 22 )* 


(m l2 n 22 - n 12 m 22 ) 
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four terms must vanish. But this can only happen for general Ui >k if each of the four 
terms in curly brackets is separately zero. This gives conditions on the entries of the 
matrices M and N for the problem to be self-adjoint. We will leave the discussion of 
these to a subsequent paper j3j, although let us note that in the case where M and N 
are real, these reduce to det M = det N. 

Therefore, for general k, we have using (|35Jl . 

k («i,o «i,fe) = f 1 det(M + NHxjkiX)), 

\n 12 m 2 2 - m 12 n 2 2\ 



which yields (I51)l . 

If, for the particular boundary conditions under consideration, 777i 2 77 22 — 77i 2 tt7 22 = 0, 
then ()49)l should be expressed not as but in an alternative way that allows the 
above derivation to be repeated. For example, if m n n 2 i — n\\m%\ 7^ 0, 

( mu nn \ (ui,k(0)\ = ( \ _ / 77112 ni2 \ fvi, k (0)\ 

^ m 21 7221 I \ui ik (l)J \det(M + NHi ik (l))J { m 22 n 22 J {vl lk (l)) , 

and we can solve for 7/1^(0) and ui, k (l). Proceeding in this way we obtain the following 
results: 

A ±fl\/T 1 ATTT f,\\ mun 2 l - 777217711 2 

det(M + M lit (l)) = - fc (m,oK,fc), "7117721 - 777217711 ^ 0, 

^ii^i,ol-U + ^n^i,o(,uj 
,,,,, Arw n » mnm 2 2 - rn 12 m 2 i 2 

det(M + NHi k (l)) = — — — — fc (7/10 Ui,fc), mnm 22 - mi 2 7772i t 0, 

mi 2 vi,o(0)* + m 1 i«i, (0)* 

A 1 TWIT 7771277,21 - 77222^11 ;2/ . . , 

det(M + iVf/ l fc (l)) = — — — fc (TJi o TJi fc), 777 i 2 77 2 l - 777 2 2 77n 7^ 0, 

mi 2 fi,o(l) - n u ui fi {0)* 



det(M + NHi k (l)) = , n22Ul1 k 2 (ui fi \u hk ) , ni 2 n 2 i - n 22 nu ^ 0, 

ni2«i,o(l) +nnWi,o(l) 



J 4-/ A/T 1 twit 777117722 - m 2 l77i2 2/ . . 

det(M + NHi k (l)) = — ——k{ui fi \ui, k ), mnn 2 2 - m 2 in 12 t °- 

ni 2 Vl,o(0)* - 777nMi )0 (l)* 

Together with (|5Tj) these equations cover all possible cases. 



B Explicit calculations for the two-component prob- 
lem 

In this appendix we will give details of the calculation required to calculate the func- 
tional determinant of the problem defined by (|58jl on the interval [—1/2,1/2] with 
boundary conditions given by the matrices M = — diag |e* M ', e~* M ', e ipl , e~ lfil ^ and 
N = I 4 . 

The evaluation of the functional determinant det' Li is based on the evaluation of 
suitably normalised solutions. In particular we need to find Ui ;k (x) = (ui )kj i(x),Ui ik; 2(x), 
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vi,k,i(x), Ui,jfc,2(a?)) such that all but one boundary condition are satisfied. We choose a 
normalisation such that 



(M + NH 1>k (l/2)) 



( u 1A1 (-l/2) \ 

«l,fe,2(~V 2 ) 
^1,mH/ 2 ) 

V V1.fc.2H/ 2 ) / 







V det(M + NH l>k {l/2)) J 



(B.l) 



This should be compared with the analogous choices (J3l?j) and (|41|) in section For 
the twisted boundary condition under consideration, this implies 



det(M + NH ltk (l/2)) = vx, k , 2 (l/2) - e~ ild v 1A2 {-l/2\ 



fB.2) 



which is used in the main text to derive (|64J). In order to derive a solution satisfying the 
condition (|B.1J) we first note that when k is not an eigenvalue (so that, in particular, 
k ^ 0), we can invert this equation to find 



/ u w (-Z/2) \ 

Ml,fc,2(-V 2 ) 

^wH/ 2 ) 

V V1.fc.2H/ 2 ) J 



( 



{M + NH 1;k {l/2)Y 









\ 



(B-3) 



V det(M + NH lik (l/2)) J 



Given the fundamental matrix H\ )k {l/2), ()B.3|) can be used to find the initial conditions 
that have to be imposed on vl\^{x) for ()B.2j) to hold. However, what actually enters 
into the final answer is /^o, the k — > value of fi jk defined by (f63|) . Therefore we only 
need to evaluate (|B.3|) in the limit k — > 0. Although (M + NH 10 (l/2))~ 1 does not exist 
(because k = is an eigenvalue) the right-hand side does exist in the limit, because the 
determinant cancels between the inverse and the entry in the column vector leaving 
only the cofactors of the matrix which are perfectly well defined in the k — > limit. To 
calculate these cofactors we only need the entries of the fundamental matrix H\^{l/2) 
at k — 0. Thus we only need to construct the normalised fundamental matrix for k = 0. 
Defining, as in section HI Yj(x) = lim k ^ Q Hj tk (x), this means we have to construct 



YAx) 



where the 4 fundamental solutions are chosen to satisfy Y"i(— Z/2) = 7 4 . The explicit 
form of these solutions is 
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yS 2) (x) 



yS 3) (x) 



yi 4) (x) 



V* 



ifil/2 



[1 — fx 2 ) — 1 + %\iz + cosh vz — ^ sinh 



vz 



V 



^--z/i(l-/i 2 )z+%- cosh vz + ^(3 - 7 /I 2 ) sinh 1/2 



-IflX 



f 2i/j, + (1 — fi 2 )z — 2i[icoshvz + i(l — 5yU 2 ) sinhz/z e* MX ^ 



V 



—2; 
2 



—z 



- sinh z/z 

V 

- sinh z/z 



^IflX 



-2i/j, + (1 — fi 2 )z + 2i\i cosh vz + ~(1 — 5/z 2 ) sinh vz 



where v 2 = 2(1 — 3/i 2 ) and z = x + 1/2. 

Due to the condition 1/2) = I4, it is easy to see that the solution with initial 
conditions (IB.3I) reads 



Yi(x) = Ml , 0il (-//2)y( 1) (x) + MlA2 H/2)yHx) + ^, 0il (-//2)y^ (x) + v 1A2 (-l/2)y?> (x). 



,(2) 



(3), 



Therefore by finding the relevant initial conditions from (|B.3|) in the k — > limit, we 
can determine the correct normalisation for yi(x). The result is given by ()66p. 



27 



References 

[I] R. P. Feynman and A. R. Hibbs. Quantum Mechanics and Path Integrals (McGraw- 
Hill, New York, 1965). 

[2] L. S. Schulman. Techniques and Applications of Path Integration (Wiley- 
Interscience, New York, 1981). 

[3] R. Rajaraman. Solitons and Instantons. (North-Holland, Amsterdam, 1982). 

[4] K. Kirsten and A. J. McKane. In preparation. 

[5] A. J. McKane and M. B. Tarlie. J. Phys. A 28, 6931 (1995). 

[6] I. M. Gel'fand and A. M. Yaglom. J. Math. Phys. 1, 48 (1960). 

[7] See, for example, E. Brezin, J-C. Le Guillou and J. Zinn- Justin. Phys. Rev. D15, 
1544 (1977). 

[8] Aspects of Symmetry: Selected lectures of Sidney Coleman. (CUP, Cambridge, 
1985). 

[9] S. Levit and U. Smilansky. Proc. Am. Math. Soc. 65, 299 (1977). 
[10] T. Dreyfuss and H. Dym. Duke Math. J. 45, 15 (1978). 

[II] R. Forman. Invent. Math. 88, 447 (1987). 

[12] R. Forman. Commun. Math. Phys. 147, 485 (1992). 

[13] D. Burghelea, L. Friedlander and T. Kappeler. Commun. Math. Phys. 138, 1 
(1991). 

[14] D. Burghelea, L. Friedlander and T. Kappeler. Integr. Equ. Oper. Theory 16, 496 
(1993). 

[15] G. Carron. Am. J. Math. 124, 307 (2002). 

[16] D. Burghelea, L. Friedlander and T. Kappeler. Proc. Amer. Math. Soc. 123, 3027 
(1995). 

[17] R. E. G. Saravi, M. A. Muschietti and J. E. Solomin. Commun. Math. Phys. 110, 
641 (1987). 

[18] R. E. G. Saravi, G. L. Rossini and M. Fuentes. J. Phys. A 25, 6743 (1992). 

[19] O. A. Barraza, H. Falomir, R. E. G. Saravi and E. M. Santangelo. J. Math. Phys. 
33, 2046 (1992). 

[20] O. A. Barraza. Commun. Math. Phys. 163, 395 (1994). 

28 



[21] M. Lesch. Math. Nachr. 194, 139 (1998). 

[22] M. Lesch and J. Tolksdorf. Commun. Math. Phys. 193, 643 (1998). 

[23] H. A. Falomir, R. E. G. Saravi, M. A. Muschietti, E. M. Santangelo and J. E. 
Solomin. Bull. Sci. Math. 123, 233 (1999). 

[24] H. Kleinert and A. Chervyakov. Phys. Lett. A 245, 345 (1998). 

[25] H. Kleinert and A. Chervyakov. J. Math. Phys. 40, 6044 (1999). 

[26] M. Bordag, E. Elizalde and K. Kirsten. J. Math. Phys. 37, 895 (1996). 

[27] M. Bordag, K. Kirsten and J. S. Dowker. Commun. Math. Phys. 182, 371 (1996). 

[28] K. Kirsten. Spectral Functions in Mathematics and Physics. Chapman & 
Hall/CRC, Boca Raton, 2001. 

[29] R. Courant and D. Hilbert. Methods of Mathematical Physics. Interscience Pub- 
lishers, Inc., New York, 1953. Volume 1, Chapter V, Section 11. 

[30] L. D. Landau and E. M. Lifshitz. Quantum Mechanics. Pergamon Press, Oxford, 
1977. Third edition. 

[31] M. Abramowitz and I. A. Stegun (eds). Handbook of Mathematical Functions. 
Dover, New York, 1965. 

[32] M. B. Tarlie, E. Shimshoni and P. M. Goldbart. Phys. Rev. B 49, 494 (1994). 

[33] M. B. Tarlie. Nonequilibrium Properties of Mesoscopic Superconducting Rings. 
PhD Thesis, University of Illinois at Urbana-Champaign, 1995. 



29 



